Nonequilibrium Steady State of Photoexcited Correlated Electrons 

in the Presence of Dissipation 



o 
o 

(N 
< 



I 

O 

o 



> 

m 
m 

cn 
o 

On 
O 



X 



Naoto Tsuji, Takashi Oka, and Hideo Aoki 

Department of Physics, University of Tokyo, Hongo, Tokyo 113-0033, Japan 

(Dated: August 1, 2009) 

We present a framework to determine nonequilibrium steady states in strongly correlated electron 
systems in the presence of dissipation. This is demonstrated for a correlated electron (Falicov- 
Kimball) model attached to a heat bath and irradiated by an intense pump light, for which an 
exact solution is obtained with the Floquet method combined with the nonequilibrium dynamical 
mean-field theory. On top of a Drude-like peak indicative of photometallization as observed in 
recent pump-probe experiments, new nonequilibrium phenomena are predicted to emerge, where 
the optical conductivity exhibits dip and kink structures around the frequency of the pump light, a 
midgap absorption arising from photoinduced Floquet subbands, and a negative attenuation (gain) 
due to a population inversion. 
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Introduction. — Recent pump-probe experiments have 
unveiled not only ultrafast dynamics of correlated elec- 
tron systems, but also a possibility to control their phys- 
ical properties and even trigger "phase transitions" with 
irradiation of pump lights. A typical example is photo- 
induced insulator-metal transition P, [J where the 
time-resolved optical conductivity [uiv), v: the photon 
energy of the probe light] exhibits a Drude-like peak in a 
low-energy region that indicates photocarriers drive the 
system into a metal. This "photodoping" opens up a 
new frontier for controlling the phase of correlated elec- 
tron materials as an alternative to chemical doping. 

Experimentally it is known for strongly correlated elec- 
tron systems that an excited state relaxes very fast (< 
1 ps) after the pumping is turned off [s], which implies 
that the system is subject to strong dissipation. Thus we 
expect that during irradiation of the pump light the bal- 
ance between pumping and relaxation is rapidly achieved, 
and a nonequilibrium steady state (NESS) emerges. To 
identify NESS is indeed a long-standing issue in nonequi- 
librium quantum statistical mechanics. For many-body 
electron problems, NESS has been mainly studied for 
quantum dot systems (e.g., Ref. 0|), while we are still 
some way from a firm understanding of NESS for bulk 
systems that might exhibit phase transitions. On the 
other hand, the relevant dissipation mechanism is still not 
clear. Time evolution of isolated correlated electron sys- 
tems has been studied for photoexcited one-dimensional 
systems [5]] or for the nonequilibrium Falicov-Kimball 
(FK) model These studies, however, do not elabo- 

rate the role of extrinsic dissipation that can be relevant 
to the long-time steady-state behavior 0]. 

Motivated by these, we pose two issues here: (i) can 
we present a general argument for NESS determined by 
dissipation in a photoirradiated electron system, and (ii) 
are there novel phenomena that emerge specifically in 
nonequilibrium? For the former we have performed an 
exact analysis of NESS for the FK model in the frame- 



work of the Floquet technique [^, as combined with 
the dynamical mean-field theory (DMFT) [llj. For the 
latter, our new finding in the observable cr(i^) character- 
istic to NESS is dip and kink structures around v VL 
(where VL is the frequency of the pump light), a midgap 
peak arising from photoinduced Floquet subbands, and a 
negative peak (gain) arising from a population inversion. 
These predictions offer experimental opportunities to re- 
alize and control the novel state with photoirradiation. 

Model. — In order to include dissipation, let us consider 
a system coupled to a heat bath with the total Hamilto- 
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As a solvable model for the heat bath, we take the 
"Biittiker probe" reservoir , 
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i?bath = E E '^P^P^^P' 
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where cj (&Jp) creates an electron (bath's fermionic de- 
grees of freedom), Sp is the kinetic energy, and Vp is the 
coupling to the mode p of the bath. The thermal bath 
is in equilibrium with temperature T. Its chemical po- 
tential is determined so that the current does not flow 
between the bath and the system. In H„iix we have in- 
cluded a contour term [the second term on the right hand 
side of Eq. ([2])] that cancels the potential shift due to the 
coupling to the bath. Let us now define a hybridization 
function T{uj) ~ J^p TrVpS{uj — Sp), and the system dissi- 
pates with a damping rate F. For simplicity, we omit the 
Lu dependence of F. 

When one starts driving the system with a continuous 
pump light at time t — 0, a transient dynamics from the 
initial equilibrium state to an excited state should first 
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occur, but then the system will relax to an NESS within 
t ~ r~^. Here we assume that (i) NESS exists, and (ii) 
NESS should be independent of the initial condition and 
correlations between the initial state and NESS, since we 
expect that those effects will be wiped out in the presence 
of dissipation (while in the absence of dissipation this is 
shown to not necessarily be the case fis'l). 

Based on these assumptions, NESS is determined as 
follows. The absence of the initial correlations allows us 
to use the Keldysh formalism. Since the pump light is 
an ac electric field periodic in time, NESS is also time- 
periodic, so that we can employ the Floquet method. 
After integrating out the bath degrees of freedom, we 
obtain the Dyson equation in a Floquet matrix form 
denoted by hats, 
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where qR^a,k (^Qm,Ao,Ko^ respectively fuU (noninter- 
acting) retarded, advanced, and Keldysh Green's func- 
tions, Yi^'^'^ are respective self-energies, 1 the identity 
matrix, and {F)mn{^) = (5mn tanh[(a; + nr2)/2T]. In 

Eq. Q, {G^° -^)mn{uj) = {u + nn + II + ir])Srnn-{ik)rnn 

with n the chemical potential, rj a positive infinitesimal, 
and ik the Floquet matrix for the noninteracting Hamil- 
tonian 



14l | . The Keldysh component (G" 
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Eq. Q, 

on the other hand, vanishes since it is proportional to irj 
while there is a nonzero dissipation term 2irF. This 
means that, although (G"^)^" contains information on 
the initial condition before the ac field is applied as well 
as on the way in which the field is turned on, it is again 
wiped out due to dissipation. NESS is thus determined 
without ambiguity by the Dyson equation ([4]). We note 
that if we stand on a phenomenological point of view, 
we are not restricted to the specific model of the reser- 
voir (21) and ([3]). Alternatively, we can start from Eq. 
interpreting F"^ as the relaxation time for a relevant dis- 
sipation mechanism such as phonons or spins. 

As one of the simplest models of correlated electron 
systems we take the FK model [3, [3 , 
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irradiated with a uniform pump light with the vector 
potential -A(t), where is the energy dispersion of a 
single-band system, U is the interaction strength, and // 
creates a localized electron. This model is known [l^ to 
undergo a metal-insulator transition as we change U. We 
adopt the model since (i) it has a simple optical excitation 
spectrum with a single charge transfer (CT) peak around 
I' = ricT ~ U (Fig. m the black curves), (ii) charge and 
spin degrees of freedom decouple from the outset so that 



we can concentrate on charge dynamics [we drop spin 
indices in Eq. ([5])] , and (iii) it can be solved exactly within 
DMFT 11 1 even out of equihbrium 0. 

Method. — To treat periodically driven correlated elec- 
trons, we can make use of the Floquet method combined 
with DMFT Il4| (an equivalent numerical technique was 
used in Ref. llTl ). Here we adopt this method for the 
dissipative case including the Keldysh component using 
Eq. (HI). We consider the hypercubic lattice with et ~ 
— t* Yl^^i cos ki/Vd in d (— > oo in DMFT) dimensions. 
For simplicity, we assume that A{t) = A{t){l, 1, . . . , 1) 
with A{t) — — E'sinfit/ri (where E is the amplitude 
of the pump hght). Then {ek)mn = ekJm-n(eE/^l) for 
m — n even, or i VkJm-n{eE/n) for m — n odd [1J| with 
— i* X^iLi sin and Jn(z) the nth order Bessel 

function. We take t* as the unit of energy. The integra- 
tion over k may be done with the joint density of states 
Pi^.v) = Ek ^(e - ^k)S{v - vk) = e~'"--"/7T Q. 

The optical conductivity is defined in such a way that 
the current change due to an infinitesimal probe light 
SEe^"^^^ (which we assume to be parallel to the pump 
light) be 5j{t) = a{v)5£e-'^* + ^V^^' ' ' )'5'^^e-''('^+"^^)*, 
and is given by Kubo-like formula [l^ 



(^H dt \ - ek-eA(i)G^it, 
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[jk{i+t/2),jk{t^t/2)])\, (6) 



with (Tq = e^/d, r = 27r/ri, G^ the lesser Green's func- 
tion, [ , ] the commutator, (...) the statistical average, 
and jk{t) = Vk-eA(t) c\{t)ck{t). In the limit i? ^ 0, we 
reproduce the equilibrium linear-response theory. 

In calculating a{v), we divide Eq. ([6]) into two parts: 
the bubble diagram and the vertex correction [see the 
inset of Fig. [T] (b)]. The former for Re (j{v) is written as 
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duj TT[vkAk{uJ + v)vkNk{uj) 



VkAk{ijj)vkNk{i^ + v)], 



(7) 



where Ak{uj) = i[G^{uj) - G^(cc;)]/2^, Nk{uj) = 
— iGfe(a;)/27r, the trace runs over the Floquet indices, 
and the bare vertex function {vk)mn = VkJm-nieE /fl) 
for m — n even, or —i CkJ-m-ni^E/ft) for m — n odd. 

The vertex correction to cr{y)^ on the other hand, ex- 
actly vanishes in equilibrium within DMFT due to the 
odd parity of Vk- However, this no longer holds for 
nonequilibrium cases, since the bare vertex Wfc_eA(t) is 
generally not parity odd. We note that the vertex cor- 
rection strongly depends on the relative polarization of 
the pump and probe lights. A virtue of the Floquet for- 
malism is that the vertex correction can be calculated ex- 
actly with the self-consistent equations for dressed vertex 
functions [l9| . In practical calculations, we introduce a 
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FIG. 1: (color online) Real part of the optical conductivity 
cr(i^) (solid lines) in units of ao for (a) Q — 1.8, (b) 2.7 and (c) 
3.3. The dashed lines illustrate (for specific values of E) the 
results without the vertex correction. The arrows indicate the 
frequency fl of the pump light. Inset in (b) depicts diagrams 
of the bubble and the vertex correction. 



cutoff for the Floquet matrix size (typically 7), for which 
we have checked the error is small enough. In the follow- 
ing we put r = 0.05, T = 0.05, and U = 3 (so the system 
is insulating), and set the system at half- filling for both 
c and / electrons. 

Numerical results. — We show the results for (t(j/) in 
Fig.[l]for three cases: (a) the frequency of the pump light 
n = 1.8 < - 17 = 3, (b) 2.7 < ncT, and (c) 3.3 > 
r^CT- When the pump light is absent {E — 0, the black 
curves in Fig.[T]), an optical gap is clearly seen in the low- 
energy region. As we increase the amplitude E, we find 
that CT peak ^ ^Ict) collapses in all the three cases, 
which is due to the bleaching effect. On the other hand, 
the low-energy behavior of cr(j^) is dramatically different 
for the three cases: In Fig.[Hb), a positive peak appears 
around ~ 0, which implies that the system is driven 
into a metallic state. Strikingly, in Fig. [TJc) we find a 
negative weight (hatched in Fig. [1]), which suggests that 
the system has some energy gain. In the case of Fig. [Tfa) 
where one-photon process is forbidden due to 17 < ficT ^ 
U, we can notice a midgap absorption around ~ 1.2 
with a positive or negative weight in ^ ~ 0. 

We also find that there are new features, besides the 
low-energy behavior, appearing around the pump fre- 
quency V ^ in all the three cases: a kink in Fig. [T]^a) 
and dips in Figs.[IJb) and[IJc). These features are unique 




FIG. 2; Left; Time-averaged density of states A{uj) (solid 
curves) and time-averaged density of occupied states N{u}) 
(shaded regions) for (a) E = 0, (b),(c) n = 1.8, (d),(e) 
Q = 2.7, and (f),(g) = 3.3 with several E. Right: The 
corresponding effective distributions /cff (lj) = N{uj)/A{uj). 



to NESS, and are not seen in equilibrium. To reveal the 
origins of the newly found features, let us first clarify 
the effect of the vertex correction by comparing the re- 
sults with and without the correction (the solid lines and 
dashed lines in Fig. [U respectively) . We find that the cor- 
rection can contribute to cr(z^) quite significantly around 
1/ ^ ^l, creating the kink [Fig. [IJa)] and dip [Figs. [Ub) 
and [IJc)] structures. Since there is no such correction 
in equilibrium, these structures are a genuine quantum 
many-body effect unique in nonequilibrium. We empha- 
size that the effect is distinct from the so-called "hole 
burning" effects observed in inhomogeneous systems [2^ . 

Other features can be well captured by the bubble di- 
agram ([7]). To analyze them, we calculate time-averaged 
local density of states (DOS) A{uj+nil.) = J2ki^k)nn{^), 
time-averaged density of occupied states N{uj + nfl) = 
(Fig- 121 left-hand panels), and an effective 
distribution fcff{oj) = N{uj)/A{uj) (right-hand panels). 

In the absence of the pump light [E = 0, Fig. ^a.)] 
we can see an insulating state equilibrated with the 
Fermi distribution. As the pump field is turned on, 
DOS changes in different manners according to which 
regime 17 belongs to: For fl ^ ^Iqt DOS hardly changes 
[Figs. [2l^d)-[2l^g)], and the energy gap remains [Figs. [2fe) 
and[2I^g)] even when the CT peak in almost disap- 
pears [Figs.^b) and^c)]. For il < flcT in Figs, [^b) 
and (He), by contrast, the DOS is modulated to have 
photoinduced midgap states, which we can assign to the 
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FIG. 3: Crossover from an insulator to a photoinduced metal: 
the weight D (open circles, joined as a guide to the eye) and 
width 7 (filled circles) against E for — 2.7. Inset is a blowup 
with a log-log plot, where the line is a fit, D — 0.335-E^. 



Floquet subbands .14] (a set of replicas of the original 
band with a spacing fl). These states are responsible 
for the midgap peak in a{i') [Fig. [Ha)] corresponding to 
one-photon process with absorption or emission energy 
ficT 1.2. 

Let us turn to the distribution (Fig. [21 right-hand pan- 
els). It is greatly modified by the pump light, where 
electrons occupying the lower band are pushed to the 
upper band absorbing the photon energy nf2, and be- 
come photocarriers. What is newly found here is that the 
nonequilibrium distribution /oft(i^) becomes highly non- 
monotonic with characteristic periodic structures (the 
period ^ Q) as clearly displayed in the right-hand panels 
of Fig. [21 Since we cannot fit /eff(w) to the Fermi dis- 
tribution with an elevated temperature, this is again a 
genuinely nonequilibrium effect, where heating picture is 
not applicable. 

The low-energy behavior of cr(i^) can also be explained 
with the distribution. For fl = 2.7 < flcT, electrons 
in the upper part of the lower band tend to be excited 
into the lower part of the upper band [Figs. [2Id) and 
[2Ue)], whereas electrons in the lower part of the lower 
band tend to be excited into the upper part of the upper 
band for n = 3.3 > Qct [Figs, [l^f) smd^g)]. In the 
former distribution, electrons prefer to absorb infinites- 
imal energy rather than emit, which is why the photo- 
carriers induce a Drude-like peak despite the energy gap 
remaining in the DOS. In the latter, a population in- 
version within each band is realized in the steady state, 
giving a negative weight in a{i'). This implies an energy 
gain, i.e., a transmitted light acquires intensity stronger 
than the incident light. A negative a{i') is also observed 
for = 1.8 < QcT and sufficiently large E [Fig. [T^c)], 
which is now explained with two-photon processes, with 
absorption energy ^ 2il > ficT- 

How does the nature of the low-energy metallic peak 
in Fig. \V[h) evolve with the pump light amplitude E7 
We fit the deviation of (t(j/) from the equilibrium one. 



5cr(i/) = (t{i') — <t{i')\e=o, to a Drude-like expression 
aot* D-f/ [n{i^'^ -f 7^)] by least squares. Here the weight 
D and the width 7 are fitting parameters, and the fit- 
ting is performed for a low frequency region (v < 0.8t*, 
here), with the rms error turning out to be < 4%. The 
result in Fig. [31 shows that D starts to grow nonlinearly, 
and is saturated around E ^ 2. We do not have a phase 
transition with singularity, but we do have a crossover 
from an insulator to a metal. 7, on the other hand, is 
nearly constant against E. This behavior is similar to 
the chemically doped FK model, where a single-particle 
excitation has a finite lifetime at T = [l6|, so that 
the photoinduced state is also a non-Fcrmi-liquid metal. 
The inset in Fig. [31 indicates D oc E"^ for small E, which 
can be understood from the third- order nonlinear opti- 
cal process, where the correction to a{i')\E=o is given as 
(5ct(z/) (X vlmx^^^ i')E'^ in terms of the nonlin- 

ear optical susceptibility x^'^'' (optical Kerr effect) (20j . 

Conclusion. — We have shown how a nonequilibrium 
steady state is determined for photoexcited correlated 
electrons (in the Falikov-Kimball model) by introducing 
the fermionic bath model. This has enabled us to pre- 
dict the features in the optical conductivity that comprise 
kink and dip structures, midgap band, as well as negative 
peaks (gain) as hallmarks of the nonequilibrium. Future 
problems include elaborations on the coherence of the 
pump light, dependence on the nature of dissipation, and 
robustness of excited states against ultrafast relaxations. 
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